Oh: 



c^ 



EUROPEAN ORGANISATION FOR NUCLEAR RESEARCH (CERN) 



CERN-PH-EP/2006-xxx 
21 April 2006 



O 

o 

(N, 

^' Test of Colour Reconnection Models 

< 

(N 



using Three-Jet Events in Hadronic Z Decays 



(N : The ALEPH Collaboration 

O 

o 

O I Abstract 

Q ■ Hadronic Z decays into three jets are used to test QCD models of colour reconnection 



(CR). A sensitive quantity is the rate of gluon jets with a gap in the particle rapidity 
q"| distribution and zero jet charge. Gluon jets are identified by either energy-ordering or 

« ' by tagging two b-jets. The rates predicted by two string-based tunable CR models, one 

implemented in JETSET (the GAL model), the other in ARIADNE, are too high and 
k> i disfavoured by the data, whereas the rates from the corresponding non-CR standard 

^ '_ versions of these generators are too low. The data can be described by the GAL model 

assuming a small value for the Rq parameter in the range 0.01 — 0.02. 



Submitted to Eur. Phys. J. C 



The ALEPH Collaboration 

S. Schael 

Physikalisches Institut das RWTH-Aachcn, D-52056 Aachen, Germany 

R. Barate, R. Bruneliere, I. De Bonis, D. Decamp, C. Goy, S. Jezequel, J. -P. Lees, F. Martin, E. Merle, 
M.-N. Minard, B. Pietrzyk, B. Trocme 

Laboratoire de Physique des Particules (LAPP), II^P^-CNRS, F-74019 Annecy-le-Vieux Cedex, 

France 

S. Bravo, M.P. Casado, M. Chmeissani, J.M. Crespo, E. Fernandez, M. Fernandez-Bosnian, LI. Garrido/^ 
M. Martinez, A. Pacheco, H. Ruiz 

Institut de Fisica dAltes Energies, Universitat Autonoma de Barcelona, E-08193 Bcllaterra 
(Barcelona), Spain' 

A. Colaleo, D. Creanza, N. De Filippis, M. de Palnia, G. laselli, G. Maggi, M. Maggi, S. Nuzzo, A. Ranieri, 
G. Raso,^** F. Ruggieri, G. Selvaggi, L. Silvestris, P. Tempesta, A. Tricomi,'^ G. Zito 

Dipartimento di Fisica, INFN Sezione di Bari, F70126 Bari, Italy 

X. Huang, J. Lin, Q. Ouyang, T. Wang, Y. Xie, R. Xu, S. Xue, J. Zhang, L. Zhang, W. Zhao 

Institute of High Energy Physics, Academia Sinica, Beijing, The People's Republic of China^ 

D. Abbaneo, T. Barklow,^^ O. Buchmiiller,^® M. Cattaneo, B. Clerbaux,^^ H. Drevermann, R.W. Forty, 
M. Frank, F. Gianotti, J.B. Hansen, J. Harvey, D.E. Hutchcroft,3°, P. Janot, B. Jost, M. Kado,^ P. Mato, 
A. Moutoussi, F. Ranjard, L. Rolandi, D. Schlatter, F. Teubert, A. Valassi, L Videau 

European Laboratory for Particle Physics (CERN), CH-1211 Geneva 23, Switzerland 

F. Badaud, S. Dessagne, A. Falvard,^" D. FayoUe, P. Gay, J. Jousset, B. Michel, S. Monteil, D. Pallin, 
J.M. Pascolo, P. Ferret 

Laboratoire de Physique Gorpusculaire, Universitc Blaise Pascal, IlsPP'^-CNRS, Clermont-Ferrand, 
F-63177 Aubicre, France 

J.D. Hansen, J.R. Hansen, P.H. Hansen, A.C. Kraan, B.S. Nilsson 
Niels Bohr Institute, 2100 Copenhagen, DK-Denmarl^ 

A. Kyriakis, C. Markou, E. Siniopoulou, A. Vayaki, K. Zachariadou 

Nuclear Research Center Dcmokritos (NRCD), GR-15310 Attiki, Greece 

A. Blondel,^2 j _q Brient, F. Machefert, A. Rouge, H. Videau 

Laoratoire Leprince-Ringuet, Ecole Poly technique, IN^P^-CNRS, F-91128 Palaiseau Cedex, France 

V. CiuUi, E. Focardi, G. Parrini 

Dipartimento di Fisica, Universita di Firenze, INFN Sezione di Firenze, 1-50125 Firenze, Italy 

A. Antonelli, M. Antonelli, G. Bencivenni, F. Bossi, G. Capon, F. Cerutti, V. Chiarella, P. Laurelli, 

G. Mannocchi,'"' G.P. Murtas, L. Passalacqua 

Laboratori Nazionah deU'INFN (LNF-INFN), 1-00044 Frascati, Italy 

J. Kennedy, J.G. Lynch, P. Negus, V. O'Shea, A.S. Thompson 

Department of Physics and Astronomy, University of Glasgow, Glasgow G12 8QQ, United Kingdom^'^ 

S. Wasserbaech 

Utah Valley State College, Orem, UT 84058, U.S.A. 

R. Cavanaugh,^ S. Dhamotharan,^^ C. Geweniger, P. Hanke, V. Hepp, E.E. Kluge, A. Putzer, H. Stenzel, 
K. Tittel, M. Wunsch^^ 

Kirchhoff-Institut fiir Physik, Universitat Heidelberg, D-69120 Heidelberg, Germany^^ 

R. Beuselinck, W. Cameron, G. Davies, P.J. Dornan, M. Girone,^ N. Marinelli, J. Nowell, S.A. Rutherford, 
J.K. Sedgbeer, J.C. Thompson,^'' R. White 



Department of Physics, Imperial College, London SW7 2BZ, United Kingdom^^ 

V.M. Ghete, P. Girtlcr, E. Kncringcr, D. Kuhn, G. Rudolph 

Institut fiir Expcrimcntalphysik, Univcrsitat Innsbruck, A-6020 Innsbruck, Austria 



18 



E. Bouhova-Thacker, C.K. Bowdery, D.P. Clarke, G. Ellis, A.J. Finch, F. Foster, G. Hughes, R.W.L. Jones, 
M.R. Pearson, N.A. Robertson, M. Smizanska 

Department of Physics, University of Lancaster, Lancaster LAl 4YB, United Kingdom^'^ 

O. van der Aa, C. Delaere,^^ G.Leibenguth,"^^ V. Lemaitre^^ 

Institut de Physique Nucleaire, Departement de Physique, Universite Catholique de Louvain, 1348 
Louvain-la-Neuve, Belgium 

U. Blumenschein, F. HoUdorfer, K. Jakobs, F. Kayser, A.-S. Miiller, B. Renk, H.-G. Sander, S. Schmeling, 
H. Wachsmuth, C. Zeitnitz, T. Ziegler 

Institut fiir Physik, Universitat Mainz, D-55099 Mainz, Germany^^ 

A. Bonissent, P. Coyle, C. Curtil, A. Ealet, D. Fouchez, P. Payre, A. Tilquin 

Centre de Physique des Particules de Marseille, Univ Mediterranee, IN^P^-CNRS, F-13288 Marseille, 
France 

F. Ragusa 

Dipartimento di Fisica, Universita di Milano e INFN Sezione di Milano, 1-20133 Milano, Italy. 

A. David, H. Dictl,^^ q Ganis,^^ K. Hiittmann, G. Liitjens, W. Manner^^ H.-G. Moser, R. Settles, M. Villegas, 

G. Wolf 

Max-Planck-Institut fiir Physik, Werner-Heiscnberg-Institut, D-80805 Miinchen, Germany^^ 

J. Boucrot, O. Callot, M. Davier, L. Duflot, J.-F. Grivaz, Ph. Heussc, A. Jacholkowska,^ L. Serin, J. -J. Veillet 

Laboratoire de lAccelerateur Lineaire, Universite de Paris-Sud, IlSpP'^-CNRS, F-91898 Orsay Cedex, 
France 

P. Azzurri, G. Bagliesi, T. Boccali, L. Foa, A. Giammanco, A. Giassi, F. Ligabue, A. Messineo, F. Palla, 
G. Sanguinetti, A. Sciaba, G. Sguazzoni, P. Spagnolo, R. Tenchini, A. Venturi, P.G. Verdini 

Dipartimento di Fisica dell'Universita, INFN Sezione di Pisa, e Scuola Normale Superiorc, 1-56010 

Pisa, Italy 

O. Awunor, G.A. Blair, G. Cowan, A. Garcia-BeUido, M.G. Green, T. Medcalf,^^ A. Misiejuk, J. A. Strong, 
P. Teixeira-Dias 

Department of Physics, Royal HoUoway Sz Bedford New College, University of London, Egham, Surrey 

TW20 OEX, United Kingdom^° 

R.W. Clifltt, T.R. Edgecock, P.R. Norton, I.R. Tomalin, J.J. Ward 

Particle Physics Dept., Rutherford Appleton Laboratory, Chilton, Didcot, Oxon 0X11 OQX, United 
Kingdom^'^ 

B. Bloch-Devaux, D. Boumediene, P. Colas, B. Fabbro, E. Langon, M.-C. Lemaire, E. Locci, P. Perez, J. Rander, 

B. Tuchming, B. Vallage 

CEA, DAPNIA/ Service de Physique des Particules, CE-Saclay F-91191 Gif-sur-Yvette Cedex, 
France^'' 

A.M. Litke, G. Taylor 

Institute for Particle Physics, University of California at Santa Cruz, Santa Cruz, CA 95064, USA^^ 

C.N. Booth, S. Cartwright, F. Combley,^^ p.N. Hodgson, M. Lehto, L.F. Thompson 

Department of Physics, University of Sheffield, ShefReld S3 7RII, United Kingdom^^ 

A. Bohrcr, S. Brandt, C. Grupen, J. Hess, A. Ngac, G. Prange 

Fachbcreich Physik, Universitat Siegcn, D-57068 Siegen, Gcrmany^^ 

C. Borean, G. Giannini 

Dipartimento di Fisica, Universita di Trieste e INFN Sezione di Trieste, 1-34127 Trieste, Italy 



H. He, J. Putz, J. Rothberg 

Experimental Elementary Particle Physics, University of Washington, Seattle, WA 98195 U.S.A. 

S.R. Armstrong, K. Berkelman, K. Cranmer, D.P.S. Ferguson, Y. Gao,^'' S. Gonzalez, O.J. Hayes, H. Hu, S. Jin, 
J. Kile, P.A. McNamara HI, J. Nielsen, Y.B. Pan, J.H. von Wimmersperg-Toeller, W. Wiedenmann, J. Wu, 
Sau Lan Wu, X. Wu, G. Zobernig 

Department of Physics, University of Wisconsin, Madison, WI 53706, USA^^ 

G. Dissertori 

Institute for Particle Physics, ETH Honggerberg, 8093 Ziirich, Switzerland. 



^Also at CERN, 1211 Geneva 23, Switzerland. 

2Now at Fermilab, PO Box 500, MS 352, Batavia, IL 60510, USA 

^Also at Dipartimento di Fisica di Catania and INFN Sezione di Catania, 95129 Catania, Italy. 

^Now at University of Florida, Department of Physics, Gainesville, Florida 32611-8440, USA 

^Also IFSI sezione di Torino, INAF, Italy. 

^Also at Croupe d'Astroparticules de Montpellier, Universite de Montpellier II, 34095, Montpellier, France. 

^Supported by CICYT, Spain. 

* Supported by the National Science Foundation of China. 

^Supported by the Danish Natural Science Research Council. 
^"Supported by the UK Particle Physics and Astronomy Research Council. 
"Supported by the US Department of Energy, grant DE-FG0295-ER40896. 

-'^^Now at Departement de Physique Corpusculaire, Universite de Geneve, 1211 Geneve 4, Switzerland. 
-'^^Also at Department of Physics, Tsinghua University, Beijing, The People's Republic of China. 
^^ Supported by the Leverhulme Trust. 

-'^^Permanent address: Universitat de Barcelona, 08208 Barcelona, Spain. 
^^ Supported by Bundesministerium fiir Bildung und Forschung, Germany. 
"'^^Supported by the Direction des Sciences de la Matiere, C.E.A. 
^* Supported by the Austrian Ministry for Science and Transport. 
i^Now at SAP AG, 69185 Walldorf, Germany 

^°Now at Croupe d' Astroparticules de Montpellier, Universite de Montpellier II, 34095 Montpellier, France. 
^^Now at BNP Paribas, 60325 Frankfurt am Mainz, Germany 
22 Supported by the US Department of Energy, grant DE-FG03-92ER40689. 

2^Now at Institut Inter-universitaire des hautes Energies (IIHE), CP 230, Universite Libre de Bruxelles, 1050 
Bruxelles, Belgique 

2^Now at Dipartimento di Fisica e Tecnologie Relative, Universita di Palermo, Palermo, Italy. 
2^ Deceased. 

26Now at SLAG, Stanford, CA 94309, U.S.A 
27Now at CERN, 1211 Geneva 23, Switzerland 
28 Research Fellow of the Belgium FNRS 
2^ Research Associate of the Belgium FNRS 

■^"Now at Liverpool University, Liverpool L69 7ZE, United Kingdom 

■^^ Supported by the Federal Office for Scientific, Technical and Cultural Affairs through the Interunivcrsity 
Attraction Pole P5/27 

■^2Now at Hcnryk Niewodnicznski Institute of Nuclear Physics, Polish Academy of Sciences, Cracow, Poland 



1 Introduction 

A description of the hadronisation of a inultiparton system requires specification of the colour 
connections among the partons. These can be modified by higher-order interference or non- 
perturbative effects in QCD, a phenomenon commonly called Colour Reconnection (CR) PJ. 
These effects can only be studied within the framework of specific models. The present interest 
in CR arises from the precise measurement of the W boson mass in W"'"W~ —>■ qqqq events in 
e"'"e~ collisions at LEP2 energies. In this fully hadronic final state possible CR effects among 
the decay products of the two W bosons contribute the largest systematic uncertainty to mw 
j2]. Direct information on CR in this channel is obtained from inter-jet particle flow studies 
J3I1] and/or from the W mass itself by comparing different jet algorithms [31 1^. Due to limited 
statistics only extreme models can be excluded. 

CR effects might also appear within a colour singlet system like the one produced in the 
reaction e+e" -^ qq(+gluons) on the Z resonance at LEPl where high statistics data are 
available. To search for CR effects one has to look for hadronic variables which are sensitive 
to the colour flow in an event. Following a proposal in Ref. [Zj, the OPAL collaboration has 
shown that gluon jets with a rapidity gap and zero jet charge, identified in events of the type 
Z — * 3 jets, provide a sensitive means to search for CR effects, and concluded that two string- 
based CR models are disfavoured by their data [Hlini- The L3 collaboration, employing angular 
gaps in the inter-jet regions of symmetric three-jet events, arrived at a similar conclusion [TH] . 
By comparing quark with gluon jets, the DELPHI collaboration found that also the standard 
colour string model (without CR) cannot adequately describe the fraction of neutral gluon jets 
with a rapidity gap [TT]. They suggest a contribution either from a colour octet neutralisation 
mechanism ^2j or, alternatively, from colour reconnection. 

In the present paper the variables proposed in |H] are used both to test CR models and to 
investigate the discrepancy reported in ^^. The data were collected by the ALEPH detector 
at LEPl. The data are compared to QCD Monte Carlo model calculations with and without 
implementation of CR and with parameters tuned to global properties of hadronic Z decays. 
Gluon jets from three-jet events are identified by either energy-ordering or by anti-b tagging. 
The rate of neutral jets is studied as a function of the rapidity gap size. Quark jets from the 
same events are used for purposes of comparison. The influence of Bose-Einstein correlations 
is also investigated. 

2 QCD models for colour reconnection 

Some basic properties of the non-perturbative CR models to be discussed in this paper are 
listed in Table HJ The first detailed theoretical study of CR was carried out by Sjostrand and 
Khoze (SK) ^^1 in the context of a possible cross-talk among the hadronically decaying W 
bosons in W pair production and its effect on the W mass measurement. The perturbative 
QCD contribution from one-gluon exchange was found to be neghgible, but there could be a 
sizeable non-perturbative contribution. Using the leading-log approximation (LLA), each qq 
system is evolved into a parton shower which, in the large- A'^c limit (where Nq is the number 
of colours), determines a sequence of colour-connected partons (quarks and gluons) which is 
used to draw the colour string. The authors developed a non-perturbative reconnection model, 
implemented ^1] in the PYTHIA generator version 6, based on the space-time overlap of the 
colour strings. The model variants SKI and SKII may be considered as two extreme descriptions 
of the colour topology. In SKI, strings are assumed to be extended flux tubes. The probability 



Table 1: Properties of colour reconnection models. 



model 


criterion of 
reconnection 


free 
parameter 


value 


effect in 
Z^qq 


SKI 


space-time overlap 
of flux tubes 


h 


0.65 


not 
implemented 


SKII 


crossing of 
vortex lines 


- 


- 


No 


ARIADNE 
ARl 


reduce total 
string length A 


-* rcco 


1/9 


Yes 


GAL 

(Rathsman) 


reduce total 
string area A 


Ro 


0.1 


Yes 


HERWIG 


reduce cluster size 
in space-time 


-t^reco 


1/9 


Yes 



to reconnect two strings is related to the overlap integral / and is given by 



Pr, 



1 — exp(— fci ■ /) 



where ki is a free parameter. In the SKII model the information is contained in thin vortex 
lines. Reconnection is assumed to take place with unit probability if they cross each other. 
The fact that the fraction of reconnected events is predicted in model SKII can be used to fix 
the parameter kj such that this fraction is the same in the two models [ki = 0.65). In the SKI 
model, CR effects would in principle also be possible within colour singlet (CS) systems like W 
or Z decays, but they have not been implemented and are thus not testable. It is remarkable 
that the SKII model predicts that such reconnections should not occur at all within singlet 
systems since the partons emerge from a single vertex. 

The ARIADNE generator, version 4, which is based on the colour dipole model, provides 
options for colour reconnection J15j. Considering all pairs of non-adjacent dipoles, reconnections 
take place with probability l/N^ if the total string length A decreases. The A measure is defined 
from the four-momenta of n colour-ordered partons: 



A 



n-l 



2=1 



log((Pi+Pi+i)V"^o) 



with ruQ = 1 GeV. In the program, Nq may be considered a free parameter, normally set to 9. 
The option ARl used in this paper enables reconnections only within colour singlet systems, 
whereas AR2 and AR3 are foreseen for reconnections also between CS systems like W pairs. 

The model due to Rathsman |l6j is based on the so-called generalized area law (GAL). The 
total area of a string is the sum of the areas of its pieces: 



n-l 



^ = ^{{Pi +Pi+lf - {^i +"^i+l)^) 
1=1 

At the end of the parton shower, pairs of string pieces are allowed to reconnect with probability 



Ro{l - exp{-bAA)), where AA = max(0, A - A^ 



(1) 



Aj-eco is the area after a colour rearrangement and the positive parameter b is one of the two 
parameters of the Lund symmetric fragmentation function. The phenomenological parameter 



Rq should be of order 1/N^. Its value has been determined by the author to be 0.10 from a 
comparison of the model to HERA data on the diffractive structure function of the proton. 
The Rathsman program is interfaced with PYTHIA(JETSET) version 5.7 [T7] . 

The HERWIG generator, version 6.1 18] or higher, offers a quite different concept for colour 
reconnection based on the space-time structure of an event at the partonic level. The relevant 
quantity is the cluster size defined as the Lorentz invariant space-time distance between the 
calculated production points of the quark and the antiquark forming a cluster. A reconnection 
among pairs of non-adjacent clusters is performed with probability Preco = 1/9 if the sum of the 
squared cluster sizes is lowered. On average, colour reconnection in this model leads to higher 
cluster masses. 

3 QCD model tuning 

Multi-hadronic final states are described by QCD Monte Carlo generators which contain a 
number of free parameters. These have to be determined from fits to data. This has been done 
extensively using event-shape, charged particle momentum and identified hadron momentum 
distributions measured by ALEPH in inclusive Z decays in order to check the overall description 
of the data and to obtain optimal values for the free parameters JHl- This section describes the 
tuning of the QCD generators when including colour reconnection. One of the effects of CR is 
that the average particle multiplicity changes slightly. For example, the average charged particle 
multiplicity, {rich), changes by -2%, -1% and +1% for the GAL, ARl and HW-CR models, 
respectively, as compared to the non-reconnected versions without changing any fragmentation 
parameters. Therefore the most important fragmentation parameters have been re-tuned for 
the colour reconnected versions of JETSET(GAL), ARIADNE and HERWIG. The tuning of 
HERWIG is described in more detail in [20j. The best fit values of the parameters are given in 
Tables |2}E] and the x^ values are summarized in Table El The results from the standard, non- 
reconnected versions are also included in the tables. For JETSET and ARIADNE, the string 
fragmentation parameters controlling the spin, flavour and type of the produced hadrons and 
the heavy quark fragmentation parameters ec, e^ are taken from Table 8 of [19j. The tuning of 
JETSET including a simulation of Bose- Einstein correlations (model BE32) is described in |2T] . 
Also given in Tables I2H11 are the predicted charged particle multiplicities in hadronic events to 
be compared to the measured value 20.91 ± 0.22 [T^ . 

The multidimensional fitting method is described in ^^1- The set of distributions used in 
the fits is listed in Table El The experimental distributions of these variables, corrected for 
detector and ISR effects, are those presented in pi| based on data collected in 1992. They are 
combined with data from the Z peak running in 1993 in order to increase the statistics to about 
1 million hadronic events. The systematic errors are those evaluated for the 1992 sample. Since 
the total errors are dominated by systematics and since the systematic uncertainties are only 
rough estimates, the x^ values given in Table El should only be considered as a relative measure 
of the fit quality. In performing the fits of JETSET(+GAL) and ARIADNE, certain regions 
are excluded where these models clearly deviate from the data (the tails of the out-of-plane 
quantities A > 0.06, Tminor > 0.20, Pt.out > 0.7 GeV and the very low Xp < 0.014). For the 
HERWIG fits the momentum spectra of non-strange and strange mesons as well as baryons are 
included in the set of distributions. 

The large overall x^ values indicate that, despite parameter tuning, it is very difficult for 
the models to reproduce the data with a 1 to 2% accuracy. The JETSET and ARIADNE 
models underestimate the Pt,out tail by up to 30% with deviations typically around 10 a. 
This distribution is much better described by the HERWIG model. On the other hand, the 



Table 2: Tuned JETSET 7.4 parameters without and with inclusion of CR. Azimuthal isotropy in 
the parton shower is assumed (MSTJ(46)=0), since this option gives a better fit to the data. Any 
parameter given without error is fixed during the fit procedure. The fraction of reconnected events is 
given in the last line. 



parameter MC name 


JETSET 


+GAL 


+GAL 


Ro PARP(188) 


- 


0.039 ±0.011 


0.10 


A (GeV) PARJ(81) 
Qo (GeV) PAR J (82) 
a (GeV) PARJ(21) 
a PARJ(41) 
b (GeV) -2 PAR J (42) 


0.312 ±0.004 
1.50 ±0.07 
0.365 ±0.009 
0.40 
0.900 ±0.018 


0.307 ±0.006 
1.57 ±0.08 
0.364 ±0.009 
0.40 
0.815 ±0.026 


0.306 ±0.006 
1.79 ±0.04 
0.362 ±0.009 
0.40 
0.724 ±0.014 


(^ch) 

f(reco) 


20.64 


20.65 
0.10 


20.67 
0.18 



Table 3: Tuned ARIADNE 4.08 parameters without and with inclusion of CR. Any parameter given 
without error is fixed during the fit procedure. The fraction of reconnected events is given in the last 
line. 



parameter MC name 


ARO 


ARl 


CR option MSTA(35) 
Ar2 PARA(26) 





1 
9 


A (GeV) PARA(l) 
PT,min(GeV) PARA(3) 
a (GeV) PARJ(21) 
a PARJ(41) 
b (GeV)"2 PARJ(42) 


0.229 ±0.003 
0.79 ±0.05 
0.358 ±0.009 
0.40 
0.825 ±0.024 


0.230 ±0.003 
0.79 ±0.02 
0.353 ±0.009 
0.40 
0.758 ±0.015 


(nch) 
f(reco) 


20.58 


20.61 
0.15 



Table 4: Tuned HERWIG 6.1 parameters without and with inclusion of CR. Two cluster model 
parameters, CLSMR and PSPLT, have been made flavour-dependent in order to better describe the 
measured B-meson fragmentation function. The D-wave meson multiplets are switched off. Any 
parameter given without error is fixed during the fit procedure. The fraction of reconnected events is 
given in the last line. 



parameter MC name 


HWO 


HW-CR 


Preco PRECO 
min. virtuality (GeV^) VMIN2 





1/9 
0.1 


A (GeV) QCDLAM 

gluon mass (GeV) RMASS(13) 

max. cluster mass (GeV) CLMAX 

angular smearing, dusc CLSMR(l) 

angular smearing, b CLSMR(2) 

power in cluster 

splitting, dusc PSPLT(l) 

power in cluster 

splitting, b PSPLT(2) 

decuplet baryon weight DECWT 


0.190 ±0.005 
0.77 ±0.01 
3.39 ±0.08 
0.59 ±0.03 


0.945 ±0.018 

0.33 

0.71 ±0.06 


0.187 ±0.005 
0.79 ±0.01 
3.40 ±0.08 
0.66 ±0.04 


0.886 ±0.017 

0.32 

0.70 ±0.06 


(nch) 
f(reco) 


20.96 


20.98 
0.08 



Table 5: x^ values for the event-shape and charged particle momentum distributions used in the 
global fits. The total x^ values for the fits using restricted regions are given in the last line. 





model 


JETSET 


GAL 


GAL 


ARO 


ARl 


HWO 


HW-CR 




CR par. 





0.04 


0.10 





1/9 





1/9 


distribution 


bins 
















sphericity, S 


23 


5 


7 


17 


27 


24 


107 


127 


aplanarity, A 


16 


91 


117 


132 


44 


60 


76 


93 


thrust, T 


21 


51 


21 


8 


15 


12 


378 


411 


minor, Tminor 


18 


76 


99 


130 


64 


81 


162 


200 


Feynman Xp 


46 


196 


171 


183 


239 


211 


370 


365 


Pt,out 


19 


1067 


1010 


1023 


864 


730 


110 


116 


Pt;m 


25 


107 


63 


70 


156 


73 


171 


164 


sum 


168 


1592 


1488 


1563 


1409 


1190 


1374 


1476 


with cuts 


137 


355 


304 


382 


372 


299 


- 


- 



distributions of sphericity and thrust, and of the scaled momenta, Xp, of charged particles 
are better described by JETSET or ARIADNE. HERWIG shows a significant excess (10 a) 
of low-thrust events and does not adequately describe the production of heavier particles like 
K-mesons and baryons. 

The total uncertainties of the parameters in Tables I2H11 are determined as follows. First, 
the systematic uncertainties of the fitted distributions are included in the x^ minimization and 
thus propagated into an uncertainty on the parameters. Second, the fit procedure is changed 
by varying the fit regions and, in the case of HERWIG, the set of baryons included in the fit. 
The largest changes in the fitted parameters are added quadratically to the first contributions. 

In order to test whether or not colour reconnection improves the overall description, the 
CR parameter of each model is fitted simultaneously with the other fragmentation parameters. 
The fit of the JETSET+GAL model has a x^ minimum with respect to the CR parameter 
Rq at approximately 0.04. If this parameter is fixed at the recommended value of 0.10, the 
changes in the other parameters are mainly a larger value for the parton shower cut-off Qo 
and a smaller value for the b parameter. The optimal value for the CR parameter A''^ in the 
ARl model is 8.4 ±1.1, consistent with the default value of 9. Although the changes in total 
X^ are not large, the distributions of global variables show a slight preference for the string- 
based colour reconnection models. For the HERWIG-CR model, the x^ increases continuously 
with increasing Preco parameter. However, the global variables like event-shape and particle 
momentum distributions considered here may not be specific enough for testing CR models. 

4 Effects of colour reconnection in Z ^ hadrons 



According to the string-based models ARl and GAL, colour reconnection leads to shorter 
strings and thus to less particle multiplicity, on average. The changes of particle multiplicity 
and of the particle momentum distributions can be compensated, at least partially, by re-tuning 
the important fragmentation parameters. Clearly, more specific variables are needed to test 
CR models. 

Studies of ARl and GAL at the generator level show that the fraction of events with one 
or more reconnections, f(reco), strongly rises with the number of partons and with the number 
of hadronic jets in the event. Also, a minimum of four partons in the final state are needed to 



perform a reconnection. This suggests that CR is related to the presence of gluon jets. The 
simplest configuration which is expected to be favourable for testing CR models is therefore a 
three-jet event with one of the jets being a well separated and energetic gluon jet. 

It is therefore important to check how well the gluon jet is described by the QCD models. 
The properties of gluon jets in comparison to those of quark jets have been studied in great 
detail at LEP [HI 122 EHl 123 ESj • These studies include the fragmentation function and its scale 
dependence, the distributions of particle multiplicity and rapidity and the sub-jet structure. 
In general, good agreement between the data and the QCD model predictions is found. An 
exception is the fragmentation function at large x where the predictions are systematically low 
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A very specific and rare class of three-jet events in which the gluon jet exhibits a rapidity gap 
and in which the hadronic system beyond the gap has zero electric charge, has been proposed by 
the OPAL collaboration 0, referring to theoretical ideas of Ref. [7j, as a signature for possible 
effects of colour octet neutralisation or colour reconnection. If a reconnection occurs according 
to the GAL or ARl models, a gluon jet, which in general consists of several gluons, often 
hadronizes as a closed string separated from the qq string by a rapidity gap. The signature is 
an increased rate of gluon jets exhibiting a rapidity gap and zero electric charge of the system 
beyond the gap. 

5 ALEPH detector, data and Monte Carlo samples 

A description of the detector and its performance can be found elsewhere [23 I2ZJ • Charged 
particles are detected in the central part consisting of a precision sihcon strip detector (VDET), 
a cylindrical drift chamber (ITC) and a large time projection chamber (TPC). Jets originating 
from heavy quarks, in particular b-quarks, are identified with a lifetime tagging algorithm which 
takes advantage of the 3-dimensional impact parameter resolution of charged particle tracks. 
The tracking chambers are surrounded by the electromagnetic calorimeter (ECAL) located 
inside the magnet coil, and the hadronic calorimeter (HCAL). 

The information from the tracking detectors and the calorimeters is combined in an energy- 
flow algorithm [2Zj. This algorithm provides, for each event, a list of reconstructed objects, 
classified as charged particles, photons and neutral hadrons, and called energy flow objects in 
the following. Objects reconstructed in the luminosity detectors are omitted. 

Multihadronic Z decays are preselected by requiring at least five good charged tracks 
whose energy sum exceeds 10% of the center-of-mass energy, -Ecm- Good tracks are defined 
as originating close to the interaction point (with transverse impact parameter |(io| < 2 cm and 
longitudinal impact parameter |2;o| < 10 cm), having at least four TPC hits and a polar angle 
such that |cos0| < 0.95. Residual backgrounds from r pair and 77 events are reduced to a 
negligible level by requiring the events to have at least 14 energy flow objects whose energy 
sum, i^vis, exceeds 50% of Ec^. This selection yields 3 378 000 hadronic events from the ALEPH 
data collected at the Z peak (-Ecm ~ 91.2 GeV) during the years 1992-1995. 

The analysis of the present paper relies on comparisons of data with QCD model 
calculations. Monte Carlo events were generated using tuned parameters as given in section 3. 
These events were passed through a full simulation of the ALEPH detector and were subject to 
the same reconstruction and analysis programs as the data. The numbers of hadronic Z decays 
generated for each of the models are given in Table IHl If in the following the GAL model is 
mentioned, the version with Rq = 0.10 is implied unless stated otherwise. 



Table 6: Statistics generated for each QCD Monte Carlo model. 



QCD model 


million events 


JETSET 


7.6 


JETSET+GAL, Ro = 0.04 


0.5 


JETSET+GAL, Rq = 0.10 


3.4 


JETSET+BE32 


0.5 


ARIADNE ARO 


3.4 


ARIADNE ARl 


3.4 


HERWIG 


0.5 


HERWIG CR 


0.5 



6 Analysis of energy-ordered jets 

6.1 Three-jet event selection 

The Durham (or A;t) cluster algorithm ^28^ is applied to the energy flow objects in order to 
determine the number of jets. The distance measure between any two particles is defined as 



Vik = 2imn{ElEl){l - cos9,k)/E^,^ = (2A;T/^vis)'. 

A value of 0.02 has been chosen for the resolution parameter ycut as a compromise between well 
separated jets and sufficient statistics. This value corresponds to fcx.cut = 6.5 GeV for the case 
Eyis = Ecm- This analysis deals with events which have exactly three jets. Their fraction is 23% 
in data before any further cuts. Three-jet events in which more than 95% of any jet energy is 
carried by photons are assumed to be hard bremsstrahlung off quarks and are rejected (0.7% 
of the three-jet events). To improve the particle acceptance, the event is only kept if the angle 
9j of each jet with respect to the beam direction satisfies | cos6'j| < 0.90. 

The jet energy resolution is improved by the following procedure. The jets are first ordered 
according to their observed energies, Ei > E2 > E3. As the measured jet vectors in general do 
not exactly form a plane, an average plane is constructed by taking the vector (pi xp3) + (p3 xp2) 
as the normal to this plane. Jets 1 and 2 are projected into this plane. The jet energies are 
recomputed from the relative angles assuming energy-momentum conservation for massless jets. 
The jets are then ordered again, but this time according to the computed jet energies such that 
El > E2 > E^. As further kinematic cuts the smallest of the interjet angles is required to be 
greater than 40 degrees and the smallest of the scaled jet energies, Xj = 2Ej/Ecm, to be greater 
than 0.1. These cuts only become important for resolution parameter values smaller than 0.02 
used for systematic checks. All these criteria result in 539 000 three-jet events (= iVsj). Some 
of the jet properties are listed in Table [3 

The information on the shower development of JETSET Monte Carlo events is used to 
estimate the probability, Pg, of a jet to be the gluon jet. The primary quarks from Z ^ qq 
(after termination of the parton shower) are assigned to the reconstructed jets by means of the 
smallest angle. The jet which has no quark assigned to it is then called the gluon jet. The least 
energetic jet (jet 3) has an average computed energy of 17.7 GeV. It represents a gluon jet in 
about 69% of the cases. The highest energy jet (jet 1) is dominantly a quark jet with fiavour 
composition given by the electroweak Z couplings. The selected sample comprises a wide range 
of kinematic configurations. 



Table 7: Some properties of energy-ordered jets in data. 



jet number 


(Ejet), GeV r.m.s. 


(^ch) 


P. 


1 
2 
3 


40.8 2.7 
32.7 4.5 
17.7 5.1 


8.81 
8.28 
7.02 


0.059 
0.248 
0.693 



6.2 Jet charge distributions 

The jet charge is defined here as the sum of the charges qi (in units of the elementary charge) 
of the particles which are assigned to jet j by the cluster algorithm: 



Qi = Y.1i 



n+ 



n. 



(2) 



i=l 



with n being the number of charged particles in the jet under consideration and n+(ra_) the 
number of particles with charge +1 (-1). This definition does not explicitly depend on the 
particle momenta. It has previously been used by the OPAL collaboration [23 in a study of 
quark and gluon jet charges. 

The Qj distribution is influenced by many factors. The charge of the underlying hard parton 
is expected to be reflected in the charges of the produced hadrons. Therefore the charge of a 
gluon jet should be zero on average. This is also true for quark-initiated jets because quarks 
are not distinguished from antiquarks in this analysis, which enter with equal frequency. An 
important property is the r.m.s. width of the distribution which depends on the mean particle 
multiplicity, on the jet environment and on the jet finder used. In addition the width is limited 
by local charge compensation in jet fragmentation together with charge conservation in the 
whole event. Also, Bose-Einstein correlations have an infiuence on the width. Experimental 
effects in general tend to smear and to shift the distribution. 

The measured Qs distribution of jet 3 is shown in Fig. [T] for the full data sample, together 
with the JETSET prediction. The data are rather well described by the simulation, in particular 
the fraction of neutral (Q3 = 0) jets. The data distribution is not symmetric around 0, but 
slightly shifted towards positive values, a feature which is also reproduced by the simulation 
and which can be explained by protons from nuclear interactions in the detector material. The 
same is true for jets 1 and 2. The colour reconnection model GAL predicts a higher rate of 
neutral jets in the gluon-enriched jet (jet 3) than is seen in the data. 



6.3 Rapidity gaps 

The particle distribution in jets is analysed in terms of rapidity 



In 



E + pl 

'm? + p^ ^ 



The longitudinal and transverse momentum components (pl and px) are measured with respect 
to the jet axis, which is defined as the vector sum of all energy fiow objects assigned to the 
jet. The pion mass is assumed for all charged particles except for those identified as electrons 
or muons, where irte or m^ is used. Neutral particles are assumed to be massless. They are 
included only if their energy exceeds 0.6 GeV because of the poor agreement between data and 
simulation at low energies. 



Figure 121 shows the multiphcity distribution of charged plus neutral particles, Uc+ni^y), for 
jet 3 within a central rapidity interval < y < y^ with the upper limit chosen as y^ = 1.5. 
About 7% of the jets in data have zero particles in this interval, as seen in the leftmost bin in 
the figure. Starting from y^ = 1.0, the rate of these 'gap-jets' falls nearly exponentially with 
increasing interval size. The rate is seen to be sensitive to CR effects : JETSET predicts too 
few and GAL predicts too many of these jets. The rate of jets with a rapidity gap obviously 
depends on the shape of the multiplicity distribution. The r.m.s. width of the distribution is 
slightly underestimated by the QCD models. 

A second definition of a rapidity gap, based on the maximum rapidity distance, A^/max, 
between adjacent particles in a jet has been proposed in jHj. Since the system beyond the gap 
selected in this way is found to consist dominantly of one particle only and is found to be less 
sensitive to CR effects, it is not investigated further. 

6.4 Charge distributions of jets with a rapidity gap 

The charge distribution of jets which exhibit a central rapidity gap as defined above is shown 
in Fig. 121 for jet 3. The jet charge is computed from charged particles with rapidities y > y.^. 
The fraction of particles moving backwards {y < 0) is negligibly small (of order 10~^). The 
distribution is normalized to the total number of three-jet events, A^sj, which means that the 
rate of gap-jets also enters the comparison of data with Monte Carlo models. The distribution 
in Fig. 121 is narrower than the corresponding one for the full sample (Fig. ^ due to the smaller 
average particle multiplicity. The rate of neutral jets with a rapidity gap (Fig.|21) shows an even 
higher sensitivity to CR than the rate of gap-jets alone. None of the models agrees well with 
the data : while the GAL model predicts too many neutral jets, standard JETSET predicts 
too few of them. 

The relative model-data differences of the rate r(0) of neutral jets 

. r(0)Mc - '^(O)data , .^^ A^(Qj = 0, Ay gap) 

5 = — , where r = '- — 3 

r(0)data A^3j 

are shown in Fig. |31 for jet 3 as a function of the rapidity gap size. Note that the points are 
not statistically independent : each point represents a sub-sample of the previous point (s). 
Both the colour-reconnected and the standard versions of JETSET and ARIADNE are seen to 
diverge from the data with increasing gap size. For the case y^ = 1.5, the GAL model deviates 
from the data by 17 cr (statistical). The CR effect in the ARl model is slightly weaker than that 
in the GAL model. Colour reconnection as implemented in these two models and assuming 
default values for the strength parameters, is therefore disfavoured by the data. The GAL 
model with a smaller value for the strength parameter, Rq = 0.04, still gives somewhat too 
high predictions. On the other hand, the rates predicted by standard JETSET and ARIADNE 
(without CR) are systematically low if there is a rapidity gap. The quantity r(0) has little or 
no sensitivity to the HERWIG-CR model, as can be seen from Fig. El 

The 6 values for the quark-enriched jet (jet 1) are shown in Fig. IHl for comparison. As 
expected, this jet shows very little sensitivity to CR for all models. The data are quite well 
described by JETSET. The ARIADNE model shows some deviation in the case of a gap. 
HERWIG is unable to describe the data. This is related to the observation that HERWIG's 
particle rapidity distributions are systematically lower over a wide range of rapidity and that 
the Qj distributions are narrower than the corresponding data distributions. 

The effect of Bose-Einstein correlations, as simulated with the tuned BE32 model, is to 
slightly broaden the charge distributions, leading to a downward shift of the 6 values by about 



-0.02 for the full samples of jets 1 and 3. For events with a gap in jet 3 (Fig. 0]) the shift in 6 
goes into the other direction. 

The results are less dependent on model-data discrepancies in the multiplicity distributions 
if the Qj distributions of jets with a rapidity gap are normalized to unit area for both data and 
Monte Carlo models. The quantity studied is the fraction of neutral jets 

N{Ay gap) 

where N{Ay gap) is the number of jets with a rapidity gap. The model-data differences are 
similar to those discussed above but statistically less significant. The predictions of the non- 
CR models JETSET and ARIADNE for the third jet are low also in this case, with a relative 
deviation from the data by about -6% for a rapidity gap in < y < 1.5, confirming the results 
of a similar analysis by DELPHI [1 1] . 

6.5 Corrected rate of neutral jets 

The results presented so far indicate that none of the QCD models provides a satisfactory 
description of the data for jet 3 with a rapidity gap. If the data are to be described by a colour 
reconnection model, the question then arises which value of the CR parameter is required by 
the data. For this study, the parameter Rq of the GAL model is considered, (Eq. (1)). In 
order to save computing time, the data are first corrected for the effects of the detector, of the 
reconstruction chain and the analysis cuts by means of a correction factor (C) derived from the 
JETSET simulation : 

m'^. = r(0)data ■ '-^^p^ = r(0)aata ■ C, (5) 

'^(0)MC,sim 

and similarly for /(O). The calculation including full simulation of the detector is denoted 
as 'MC,sim'. The calculation at the particle level of the event generator ('MC,gen') is done 
separately according to the following prescription. Hadronic Z events are generated without 
initial- and final-state photon radiation. The Durham cluster algorithm is applied to all charged 
and neutral stable particles to select three-jet events. The jets are ordered according to their 
actual energies. Rapidities are calculated using the actual particle masses. To define a rapidity 
gap, neutral particles with energies smaller than 0.6 GeV are omitted. 

Table 8: Corrected data values for the rates of neutral jets in jet 3 with a rapidity gap from to 1.5. 
The first error is statistical, the second systematic. 



r(0) 
/(O) 



corrected data 



0.0234 ± 0.0004 ± 0.0015 
0.437 ± 0.006 ± 0.010 



The results for the lowest energy jet (jet 3) and requiring a gap in the rapidity range 0-1.5 are 
given in Table |H1 The quoted systematic error arises from model dependence and is estimated 
by taking the largest difference of the results when correcting with JETSET, GAL(0.04) or 
ARIADNE. These numbers can then be compared to those obtained from generator level 
calculations of the JETSET-I-GAL model for different values of Rq. The complication here 
is that the CR effect is not a one-parameter problem : the fraction of reconnected events 
depends strongly on the parton shower cut-off, Qq, but only very weakly on the fragmentation 



parameter b. To ensure agreement of the average particle multiplicity with the data, at each 
Rq value the parameters A, a and b are re-tuned to global quantities, while keeping Qq fixed 
at 1.5 GeV. 

The comparison yields an optimal value for the colour reconnection parameter Ro of about 
0.02 if the quantity r(0) is used, and about 0.01 if /(O) is used. 

7 b-tag analysis 

A much higher gluon jet purity can be achieved by identifying both the quark and the antiquark 
jets in events of the type Z— >^ bbg, by utilizing the long lifetime of B-hadrons [i^Q,. Based on 
the three-dimensional impact parameter and its significance measured for each charged particle 
track the probability Pjet is computed that all charged tracks of a jet arise from the primary 
interaction vertex. The method of gluon-tagging is taken from J22I- Three-jet events are selected 
in which two of the jets exhibit significant lifetime (Pjet < 0.01), whereas the remaining jet does 
not (Pjet > 0.01). This latter jet then represents the tagged gluon jet. No energy ordering is 
applied here. 

The event selection and the computation of the jet energies are the same as already described 
in Sect. 6.1 except that all three jets are required to fall in the geometrical acceptance of the 
VDET (|cos6'j| < 0.766) and that the scaled gluon jet energies, Xg, are restricted to values 
smaller than 0.85. 

The fraction of three-jet events surviving the lifetime tag condition alone is 7.6% in data, 
to be compared to 7.4% in the JETSET simulation. About 24600 tagged three-jet events are 
selected from data. According to JETSET Monte Carlo calculations 95% of the events actually 
arise from primary b-quark pairs, 4.2% from c-quark pairs and the rest from light quark pairs. 
The tagged gluon jet energies are distributed with mean value 19.8 GeV and an r.m.s. of 7.8 
GeV. The tagged gluon jet is estimated to be the 'true' gluon jet in 97.2% of the cases on 
average. Since the gluon purity is found to drop rapidly at high energies, the scaled jet energies 
are restricted to values Xg < 0.85 as mentioned above. A sub-sample of 18700 tagged gluon 
jets are lowest energy jets. This sub-sample is therefore in common with the jet-3 sample of 
the energy-ordered analysis presented in Sect. 6. 

In order to test colour reconnection models the same analysis is performed on the tagged 
gluon jet as on the least energetic jet in Sect. 6. Due to the higher gluon purity the results 
are expected to be less dependent on the quark jet background in the sample. Figure [7| shows 
the charge distribution of the tagged gluon jets. The rate of neutral (Qg = 0) jets is well 
described by JETSET. The CR model predicts a clear excess with respect to the data. Figure 
IHl shows the multiplicity distribution of charged plus neutral particles within the central rapidity 
interval <y < 1.5. The rate of jets having zero particles in this interval (rapidity gap) shows 
sensitivity to CR. 

The charge distribution of the jets exhibiting a rapidity gap is shown in Fig. El This 
distribution is normalized to the number of tagged three-jet events and contains 1090 jets in 
data. Clearly, the CR models GAL and ARl predict too many neutral gluon jets. The relative 
model-data differences 6 (Fig. ITU]) are numerically larger than those from the energy-ordered 
analysis and reach values around 1. The results of the b-tag analysis confirm those of the 
energy-ordered analysis although with less statistical significance. At y^ = 1.5 the GAL model 
deviates from the data by about 7 a (statistical). 

The distribution of the charged track multiplicity in gluon jets with a rapidity gap is shown 
in Fig. m The data do not exhibit clear evidence for spikes at the even values n^h = 2,4 and 
6 as predicted by the CR models GAL and ARl. 



Table 9: Experimental cuts varied for systematics studies. 



cut 


standard value 


changed value 


\do , (cm) 


2.0 


0.5 


kol, (cm) 


10.0 


5.0 


Ntpc 


4 


6 


cos 9j 


0.90 


0.80 



8 Systematic checks 

The energy-ordered analysis has been repeated by varying the definition of the jets and the 
definition of the rapidity gap in the following ways : 

• The Durham jet resolution parameter, Hcnti has been varied in the range from 0.01 to 
0.03. 

• A different jet finder, the JADE algorithm, has been used with resolution parameter 
2/cut = 0.08, which leads to a similar three-jet event rate. 

• A re-assignment of particles to jets was performed on the basis of the smallest angle with 
respect to the jet axes. Afterwards the jet four-momenta are re-computed. This mainly 
affects particles at large angles to the jet axes which are of importance for defining a 
rapidity gap. 

• The rapidity gap has been defined with charged tracks only. 

• The shifted rapidity regions (0.5-1.5) and (1.0-2.0) were also used to define the gap. 

In all five cases the qualitative features of the results remain the same, although the numerical 
values of the measured quantities may change. 

Possible inadequacies of the detector simulation have been studied by varying, in data and 
simulation, some of the cuts on the charged tracks and the cut on the polar angle of the jet 
axes according to Table IHl This has been done only for the normalized jet charge distributions. 
In all cases, the model-data differences only change within their statistical errors. 

The model-data comparison of the rate of jets with a rapidity gap relies on the calibration 
of the average particle multiplicity at the detector level. In all hadronic events, the average 
charged particle multiplicity, (rich), of all generators, including HERWIG, agrees with the data 
to an accuracy below 1% which, of course, is a result of the model parameter tuning. For 
JETSET and ARIADNE this also holds for three-jet events, and in addition for each of the 
three jets individually. In contrast, HERWIG exhibits a 1.5 — 2.5 % deficit in each of the three 
jets of three-jet events, thus giving less reliable predictions. 



9 Conclusions 

Three-jet final states from e"'"e~ anniliilation into liadrons are used to test QCD models of colour 
reconnection (CR). From a sample of approximately 3.4 million hadronic Z decays, collected 
by the ALEPH detector at LEP, three-jet events are selected by the /ct cluster algorithm with 
a fixed cut-off, ?/cut = 0.02. The distributions of particle multiplicity in jets in fixed rapidity 
intervals and the distributions of the jet charge are investigated. The main analysis is based 
on jet energy ordering. According to JETSET, the least energetic jet (jet 3) represents a gluon 
jet in about 69 % of the cases. 

Gluon jets with a large gap in rapidity and with zero electric charge of the system beyond 
the gap constitute a particularly sensitive test. The rate, r(0), of these jets is measured as 
a function of the rapidity gap size and is compared to predictions of globally tuned QCD 
generators. The rate r(0) is predicted too high by the string-based CR models JETSET+GAL 
and ARIADNE ARl when using default values for the CR strength parameters. Thus these 
two models are disfavoured by the data as already observed in Ref. [Hj. On the other hand, 
the rate is too low for the standard, non-CR versions of these generators, in agreement with a 
similar analysis in Ref. JT] . The Ro parameter of the GAL model is constrained by the data 
to the range 0.01 — 0.02. The quantity r(0) turns out to be insensitive to CR as implemented 
in the HERWIG-CR model. 

In a separate analysis of the same data, two b-quark jets are positively identified using 
lifetime tagging, thus achieving a very high gluon jet purity of about 97 % for the remaining jet. 
The results confirm those of the energy-ordered analysis though with less statistical significance. 

Assuming that the physics of colour rearrangements is similar in hadronic Z decays and in 
WW decays, this result would imply that the CR models as implemented in JETSET-(-GAL 
and ARIADNE overestimate the systematic shift of the measured W boson mass. 
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Figure 1: The charge distribution of jet 3 (the lowest energy jet) compared to JETSET predictions 
without and with colour reconnection. 
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Figure 2: The distribution of the number of charged plus neutral particles in the rapidity interval 
< y < 1.5 for jet 3, compared to JETSET without and with CR. 
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Figure 3: The charge distribution of jet 3 for a rapidity gap in < y < 1.5. 
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Figure 4: The relative model-data differences 6 in the rate r(0) of neutral jets for jet 3 as a function 
of the upper limit of the rapidity gap. The points at y^ = correspond to the full sample. Statistical 
uncertainties of the data and the Monte Carlo calculations are shown as error bars. 
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Figure 5: The relative HERWIG-data differences 5 in the rate r(0) of neutral jets for jet 3 as a 
function of the upper limit of the rapidity gap. 
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Figure 6: The relative model-data differences 6 in the rate r(0) of neutral jets for jet 1 as a function 
of the upper limit of the rapidity gap. The points at y„ = correspond to the full sample. Statistical 
uncertainties of the data and the Monte Carlo calculations are shown as error bars. 
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Figure 7: The charge distribution of the anti-b tagged gluon jet compared to JETSET predictions 
without and with colour reconnection. 
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Figure 8: The distribution of the number of charged plus neutral particles in the rapidity interval 
< 2/ < 1.5 for the gluon jet, compared to JETSET without and with CR. 
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Figure 9: The charge distribution of gluon jets with a rapidity gap in < y < 1.5. 
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Figure 10: The relative model-data differences 5 in the rate r(0) of neutral jets as a function of the 
upper limit of the rapidity gap for the gluon jet sample. The points at y^ = correspond to the full 
sample. Statistical uncertainties of the data and the Monte Carlo calculations are shown as error bars. 
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Figure 11: The charged multiphcity distribution of gluon jets with a rapidity gap in < y < 1.5. 



